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SUMMARY 

This  Report  summarises  the  available  knowledge  of  the  component 
losses  in  a combustion  chamber.  The  information  given  in  this  Report  should 
enable  the  pressure  drops  through  swirlers,  primary  baffles,  cooling  systems, 
etc,,  to  be  calculated.  Most  of  the  data  were  abstracted  and  collected 
from  the  various  reports  listed  in  the  bibliography.  In  certain  cases  (e.g. 
mixing  losses)  the  information  is  incomplete  and  in  these  circumstances  the 
limited  experimental  results  available  are  supplemented  by  hypotheses  which 
require  proof,  A specimen  calculation  of  the  pressure  drop  and  air  flow 
distribution  of  a typical  chamber  is  given  in  Appendix  II.  Vhe  calculated 
and  measured  values  of  pressure  drop  (cold)  agreed  within  k per  cent. 
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1 .0  Introduction 


Effective  combustion  chamber  design  and  development  requires  a know- 
ledge of  the  air  flow  distribution  throughout  the  chamber.  Since  the  air 
flows  through  the  chamber  in  two  or  three  principal  paths,  arranged  in 
parallel,  the  loss  of  total  pressure  in  each  path  must  be  the  same.  Thus 
the  division  of  air  between  the  various  paths  will  be  determined  by  their 
relative  resistances.  This  resistance  to  flew  in  each  path  is  the  summation 
of  the  individual  component  losses.  For  example,  the  primary  cirouit  resis- 
tance comprises  the  swirler  loss,  diffusion  loss,  oembustion  loss,  eto. 

Hence  to  obtain  the  air  distribution  in  a given  ohamber  the  component  losses 
must  be  calculable  from  design  dimensions,  and  a method  of  combining  the 
oirouit  resistances  available.  Such  a method  was  developed  by  Probert  and 
Kielland'  and  subsequently  simplified2  by  dispensing  with  the  "step-by-step" 
system  of  calculation.  However,  no  comprehensive  report  on  component  losses 
has  yet  appeared  although  a note  for  disoussion  was  published^.  The  present 
Report  supplies  the  hitherto  missing  data  much  of  which  was  obtained  from 
unpublished  work  at  N.G.T.E.  In  oaseB  where  the  information  is  inoonplete 
the  available  data  are  supplemented  by  hypotheses  which  require  proof. 

In  the  Report  each  component  is  considered  in  detail  and  the  method 
of  obtaining  the  overall  loss  and  air  distribution  added  for  corrpletenesa . 
Appendix  II  gives  a specimen  calculation  for  a conventional  ohamber. 

2.0  gwirlera 


Flow  conditions  at  outlet  from  a swirler  vary  along  the  blade  span 
to  satisfy  radial  equilibrium  as  shown  in  Appendix  III.  Thus,  free  vortex 
blading  gives  a constant  axial  velocity  oomponent  while  the  whirl  velocity 
varies  inversely  as  the  radius.  Other  forms  of  blading  each  have  their  own 
particular  characteristics.  Although  true  mean  values  of  the  velocity  com- 
ponents should  be  used  for  pressure  I033  calculations,  negligible  error  is 
involved  and  the  tedium  of  obtaining  these  values  obviated,  by  using  values 
ooourring  at  the  weighted  mean  radius  (rm) . 


rm  = 


1 

/2 


(R* 


+ r. 


1 

,2 


(D 


where  the  symbols  have  the  significance  given  in  Appendix  I. 

It  is  possible  to  study  theoretically  the  efficiency  of  swirlers  in 
turning  the  air  through  a given  angle  by  considering  the  two  dimensional 
flow  of  a perfect  fluid  through  a lattioe  of  plates.  This  problem  has  been 
studied^  and  the  results  applied^  to  connect  the  angle  of  deviation  (a)  with 
the  pitch-ohord  ratio  (or)  for  various  angles  of  stagger  (3) . In  Figure  1 
angle  of  deviation  is  plotted  against  pitch-chord  ratio.  The  ourves  show 
that  for  quite  praotioal  pitch-ohord  ratios,  i.e.  0.5  < O’  < 1 .0  the  devia- 
tion angle  is  almost  identioal  with  the  stagger.  Experimental  results  agree 
with  this  finding  and  it  is  now  usual  to  employ  pitoh-ohord  ratios  of  about 
0.7  for  all  swirlers  required  to  give  a tight  swirl  (i.e.  high  values  of  a 
and  3) • Thus  for  theoretical  oaloulations  on  swirler  pressure  losses  it  is 
both  convenient  and  justifiable  to  assume  that  the  air  is  deviated  through, 
the  entire  stagger  angle  3. 
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2.1  Pressure  drop  due  to  a swirler 

By  considering  in  some  detail  the  flow  through  the  swirler  and  the 
resultant  motion  of  the  air,  an  expression  for  the  pressure  drop  oan  be 
derived. 

Consideration  is  now  given  to  the  outlet  flow  from  the  svdrler  at  the 
mean  radius  as  defined  by  Equation  ( 1 ) . 

2.11  Whirl  velooity  oanponent  dissipated  and 
oonatant  static  pressure 

Dissipation  of  the  whirl  velocity  head  is  the  most  obvious  assumption 
regarding  swirler  pressure  drops . But  an  assumption  must  then  be  made  about 
the  static  pressure  relationship  at  the  swirler  outlet  (1)  and  at  a plane  (2) 
situated  downstream  in  the  flame  tube.  A likely  assumption  is  that  the  mean 
static  pressure  difference  is  negligible. 

A mere  statement  of  the  total  pressure  loss  is  obtained  by  applying 
Berroullis  equation,  thus 


= ^2  + losa 


(2) 


with  the  further  assumption  of  constant  static  pressure  this  reduces  to 


loss  = 


(3) 


and  ainoe 


V1  = va  aec  a 


A 

Vg  = Va  (whirl  component  lost) 


% 
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2.12  Whirl  velocity  component  dissipated  and 
axial  momentum  conserved 


A more  logioal  assumption  than  constant  static  pressure  is  conserva- 
tion of  axial  momentum.  Even  this  must  have  certain  limitations  since  the 
axial  momentum  is  unevenly  distributed  aoross  the  flame  tube  diameter  and  is 
negative  along  the  axis  due  to  flow  reversal. 

The  momentum  equation  is 


Pi  a8-p2  v(^~)  (%-as)  = va2{^  -tj  (6) 


thus 


Pi 


- P2  = 


and 


$F  = 


1 + 4 


(AF\(Aa  ~ V\ 
Us/Wa  + A?) 


(7) 


2.13  Consideration  of  most  reliable  assumption 

Of  these  two  views  the  former  has  proved  to  be  the  more  reliable . 
Although  there  is  the  possibility  of  some  slight  pressure  recovery  by  virtue 
of  the  area  change  it  is  undoubtedly  local  and  is  dissipated  by  the  friotion 
in  the  ensuing  recirculation  and  general  oombustion  turbulenoe.  The  compari- 
son is  good  between  measured  losses  given  in  Reference  5 and  by  calculation 
using  Equation  (12)  whioh  is  Equation  (5)  plus  the  blade  loss.  For  typical 
values  of  Ap,  Aa  and  a the  difference  in  loss  factor  by  using  Equations 
(5)  and  (7)  rarely  exoeeds  5 per  cent,  the  former  giving  better  agreement 
with  experimental  results.  Conservation  of  axial  and  angular  momentum 
considerably  increases  the  difference  between  calculated  and  experimental 
results. 


2.14  Blade  losses 

In  the  foregoing  analysis  the  losses  are  assumed  to  originate  from 
the  resultant  flew  conditions  of  the  air  after  leaving  the  swirler  and  no 
mention  was  made  of  the  losses  in  the  swirler  itself.  These  are  due  to 
profile  and  secondary  losses  in  the  blades.  The  former  are  losses  attri- 
butable to  sldn  friction  and  separation,  the  latter  due  to  three  dimension- 
al effeots.  These  losses  are  approximately  of  the  same  magnitude  and  in  the 
oaae  of  swirlers  where  the  incidence  is  zero,  the  principal  factors  effect- 
ing the  overall  blade  loss  are  outlet  angle,  pitoh-ohord  ratio  and  blade 
passage  area.  However,  since  the  blade  loss  represents  a very  small  propor- 
tion of  the  total  swirler  loss,  an  average  figure  of  15  per  cent  of  the 
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swirler  outlet  velooity  head  is  taken  - for  the  blade  loss  for  values  of  (a) 
in  the  range  65°  < a < 85°.  This  figure  was  experimentally  determined  (see 
Reference  6)  and  is  independent  of  blade  form. 

For  smaller  values  of  (a)  and  for  increased  accuracy  where  such  vari- 
ables as  blade  height  and  thickness  are  taken  into  account,  the  following 
method  abstracted  from  Reference  7 is  used. 

This  method  is  used  for  determining  the  losses  in  turbine  nozzle 
guide  vanes  and  there  are  obvious  limitations  when  it  is  applied  to  swirlers. 
Errors  are  most  likely  to  be  associated  with  the  secondary  loss  ooeffioient. 

Hub  ratios  for  turbines  are  of  the  order  0.8  whereas  for  swirlers  they 

are  about  0.2.  Reducing  the  hub  ratio  undoubtedly  increases  the  secondary 
loss  for  turbines  and  will  presumably  affeot  swirlers  similarly,  although  to 
a greater  degree.  However,  the  deflection  angles  and  flow  accelerations  are 
higher  in  swirlers  and  the  latter  at  least  will  tend  to  reduoe  the  loss. 

These  various  effeots  are  allowed  for  (see  Section  (b)  below),  but  the  over- 
all impression  is  that  the  method  of  Reference  8 when  applied  to  swirlers 
tends  to  underestimate  the  secondary  loss.  Unfortunately  there  are  not 
sufficient  swirler  tests  for  an  independent  estimate  of  the  secondary  loss 
to  be  made. 


Conditions  are  considered  at  the  reference  radius  rm. 

Details  required  ( see  Appendix  I and  Figure  2) . 

(1)  Blade  chord  - c at  reference  radius 

(2)  Blade  pitch  - s at  reference  radius 

Blade  thickness  ~ t - at  reference  radius 

2 2 

Free  swirler  area  ~ As  = n (R  - r0  ) 

(a)  Profile  loss  coefficient 

From  Figure  3 knowing  (a)  and  the  pitoh-chord  ratio  (<r)  the  profile 
loss  ooeffioient  (Yp)  is  obtained. 


tan  a (8) 

tan  a • oos  an  ..  . - (9) 


(b)  Seoondary  loss  coefficient 

For  zero  incidence  and  assuming  a = 3 

. 1 
tan  an  = 2 

also  Ci/(s/c)  a 2 

The  seoondary  loss  for  zero  incidence 


Ys 


= K 


2 2 

[ci/(s/o)J  joos2o/oos^an 


] 


(10) 


The  factor  K is  a funotion 


and  is  plotted  in  Figure  4. 
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2 2 
Ag  = n (R  - r0  ) 

At  s An  ooe  a (An  = swirler  outlet  area) 

(o)  Total  loss  ooeffioient  (Y^) 


Yt  - Yp  + Ya  (11) 

If  the  t/s  ( thioknesa/pitoh)  ratio  differs  from  0.02  then  the  total 
loss  ooeffioient  should  be  oorreoted  by  the  multiplication  factor  given  in 
Figure  5* 


2.15  Effect  of  Reynolds  number  on  blade  losses 

The  Reynolds  number  for  a swirler  is  defined  in  the  usual  manner 
using  the  blade  chord  as  the  scalar  length  and  the  outlet  absolute  velooity, 
density  and  viscosity  at  the  mean  radius  rm.  For  all  farms  of  aerodynamic 
machine  the  loss  increases  with  decrease  of  Reynolds  number  especially  in 
the  range  Re  < ICk.  The  effect  of  Reynolds  number  on  profile  loss  may  be 
determined  approximately  from  Figure  6 which  has  relative  loss  ooeffioient 
Y 

(defined  as  ■ ■■■- — ..  ) plotted  against  Re  and  is  for  all  forms  of 

Yp  at  Re  = 2 x 1(P 

blading.  The  secondary  losses  are  as Burned  to  be  independent  of  Reynolds 
number". 


2.16  Overall  loss  ooeffioient  for  a swirler 


From  Equation  (7)  and  Section  2.1  A the  total  loss  coefficient  for  the 
swirler  in  terms  of  the  flame  tube  area  Ap. 


or  a little  more  accurately 

2 


sec  a - 1 


(12) 


=(?)  +*t}-1 


(13) 


2.1 7 Overall  loss  ooeffioient  for  various 
types  of  swirler 

Equations  (12)  and  (13)  are  quite  general  equations  for  conventional 
swirl ors  and  it  only  remains  for  one  or  two  general  observations  to  be  made 
when  these  formulae  are  applied  to  the  various  types  of  swirler. 

2.18  Constant  blade  angle  - curved  blades 

This  type  of  swirler  is  frequently  used  where  'tight'  swirls  are 
required  and  where  the  velocities  are  relatively  high.  The  blades  are  ourved 
so  that  the  upstream  edges  are  parallel  to  the  flow,  i.e,  zero  incidence. 
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Either  Equations  (12)  or  (13)  may  be  used,  to  determine  the  loss  factor  (ftp). 

2.19  Constant  blade  angle  - straight  blades 

This  type  of  swirler  is  very  easily  manufactured  and  is  representative 
of  the  swirlers  used  in  large  industrial  type  chambers,  where  the  overall 
velocities  are  low.  Sinoe  the  inoidenoe  of  the  blades  is  extremely  high,  the 
loss  factor  is  also  very  high,  although  the  flat  blades  are  extremely  effec- 
tive in  deviating  the  air  through  the  required  angle.  Scanty  evidenoe5 
suggests  that  the  blade  loss  is  approximately  doubled  compared  with  ourved 
blades  for  the  same  value  of  (a)  where  65°  < a < 85°. 

Hence  loss  coefficient  for  swirler  with  flat  plates  is  given  by 


Obviously,  the  more  accurate  oaloulation  of  Yt  i0  impossible  in  this 
instance  sinoe  the  blades  are  permanently  stalled  due  to  the  very  high 
inoidenoe . 


2.20  Varying  blade  angle  - ourved  blades 

In  view  of  the  manufacturing  difficulties  and  the  small  inorease  in 
performance  over  the  oonstant  blade  angle  type,  this  type  of  swirler  is  now 
rarely  used.  The  blades  are  usually  of  free  vortex  form  giving  maximum  whirl 
velocity  and  hence  low  pressure  at  the  oentre.  To  apply  the  loss  coefficient 
formula  it  is  neoeasary  to  ascribe  a value  to  (sec  a) . As  mentioned  in 
Section  2.0  negligible  error  is  involved  by  applying  values  occurring  at  the 
mean  radius  (r^) , As  shown  in  Appendix  III  if  the  blades  are  of  free  vortex 
form  (Va)  is  oonstant  and 


sec  a = 1 + 


£ ro 
R2  + r, 


- tan  a0 


and  ftp 


2 r„ 


R*  + rr 


tan 


1 


(15) 


For  forced  vortex  blades  (rarely  used) 


8eo  a ■ 1 + 


2 2 
R + rQ 


2(r  2 ooseo2a0-R2) 


and  ftp  ■ 


£){ 


1 + 


2 2 

R + ro 


2(r02  ooseo2  aD  “ R2 


} <1  ♦ *t> 


- 1 


(16) 
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2.21  Ported  swirler  (Figure  7) 

The  development  of  a combustion  ohamber  having  low  wall  temperature* 
resulted  in  a stabilising  baffle  embodying  this  type  of  swirler.  Assuming 
no  pressure  reoovery  and  the  static  pressure  difference  between  the  swirler 
outlet  and  flame  tube  are  negligible,  as  in  Seotion  2.13; 

loss-*  £ p (V-,2  - V22) 


From  the  velocity  triangle  of  Figure  7 


V,  = V ooseo  a 

va  = v ain  9 


loss 


1 

7 


2fl  2 
oosec  6 coseo  a 


2 

coseo 


1 


(17) 


From  Equation  (17),  as  the  semi-angle  of  the  cone  and  the  air  angle 
through  the  ports  relative  to  the  tangent  at  the  ports  increase,  the  loss 
decreases.  This  is  to  be  expected.  There  are  no  known  experimental  results 
from  which  an  allowance  for  blade  loss,  i.e.  air  friction  at  the  ports,  can 
be  made. 


2.22  Tangential  port  swirler 

This  type  of  swirler  was  last  used  on  the  early  types  cf  chamber  for 
the  W2B,  W2/500  and  W2/700  engines,  and  may  not  be  used  in  the  same  form  again. 
For  the  purpose  of  determining  the  loss  it  is  reasonable  to  assume  that  the 
velooity  head  through  the  ports  is  lost. 


2 


2.23  Vortex  type  swirler  (Figure  8) 

This  type  of  swirler  is  basically  a small  vortex  chamber  followed  by 
a throat  and  is  a comparatively  new  type . Its  ability  to  ' run  full*  gives  it 
an  advantage  over  the  conventional  swirler.  With  reference  to  Figure  8,  the 
pressure  loss  comprises  two  principal  components.  Firstly,  that  due  to 
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producing  the  whirl  velocity  at  the  throat  and  seoondly,  the  production  of 
the  axial  velooity  component. 

The  pressure  drop  between  the  tangential  entry  and  the  throat  is 
mainly  a friction  drop  and  assuming  the  vortex  decay  law 


V„  rn  = o 


Total  pressure  drop  AP  oan  be  shown  to  be 
2 2n  2n 


(19) 


by  integrating  the  equation  for  static  pressure  drop  in  vortex  flow: 
dr  ^ r 


between  r^  and  r2  and  since  the  swirl  energy  at  the  throat  is  irrecoverable 
2 2n  2n 

r2  / \r1 


"■W 


and  also  the  axial  outlet  velocity  must  be  produoed. 

Hence  total  pressure  drop 

2 2n  2n  2 

2 1 

v2  +I 


('-0} 


also  by  continuity  Vw  Ag  = Va  =»  Ag,  Vj, 


. . AP  « 


2n 


9 ;-£->-©K)i 


2 2-i 

A»  \ . 


thus 


r 2n 


(20) 


(21) 


..  ..  ..  (22) 
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For  a free  vortex  n = 1,  but  teats  on  a model  oyolone  of  approximately 
18  in.  maximum  diameter  and  6 in.  wide  have  shown  n = 0.95,  and  that  n 
deoreases  further  as  the  width  is  reduoed.  Sinoe  for  a praotioal  Bize  of 
swirler  the  effective  Reynolds  number  is  lower  and  the  flow  area/wetted  area 
ratio  is  greater,  n will  probably  be  of  the  order  0.8.  No  experimental 
results  are  available  for  confirmation  of  this  value.  The  angle  of  swirl  at 
the  throat  (w)  is  given  by 


u = tan 


-1  ^W2 


n . 

ri\  h 


2'  As 


Thus  a wide  latitude  is  allowed  in  designing  a vortex  swirler  for  a 
given  value  of  swirl  angle. 

2.3  Swirler  followed  by  a throat 

The  ccirbination  of  a swirler  followed  by  a throat  ocours  frequently 
in  chambers  containing  ceramic  liners.  This  problem  was  studied  ° and  pre- 
dicted values  for  the  pressure  loss  were  closely  substantiated  by  experi- 
mental results , The  problem  is  complicated  by  the  fact  that  heat  addition 
occurs  at  the  reference  planes  downstream  from  the  swirler  exit.  With 
reference  to  the  notational  diagram  Figure  9 


Axial  velooity  from  swirler  ^ 


outlet 


Whirl  velocity  from  swirler  Ap  f1' 


outlet 


Vp  tan  a 


The  kinetio  energy  changes  between  the  plane  of  the  swirler  outlet 
and  the  oeramio  liner  throat  are  baaed  on  the  assumptions  that  the  axial 
velooity  oomponent  increases  in  the  ratio  of  the  areas  and  the  whirl  velooity 
in  the  square  root  of  this  ratio,  making  it  a type  of  free  vortex.  This 
latter  assumption  implies  that  the  moment  of  momentum  is  oonstant  on  a stream 
surface  and  is  described  in  Reference  5. 


Ap  t" 

Axial  velooity  at  throat  = — • Vp  . — 

At  H 


Ap  /A  o' 

Whirl  velooity  at  throat  = — • /—  — • Vp  tan  a 

Ag  I At  P" 
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Thus  assuming  no  pressure  recovery  and  the  statio  pressure  difference 
to  be  negligible  between  the  throat  and  the  flame  tube  downstream 


the  blade  loss  in  the  swirler  is 


Overall  loss  factor  obtained  by  combining  Equations  (28)  and  (29)  and 
simplifying 


= 


£)(£)  I 


[iSt»2a 

I A* 


2 

sec  a 


(30) 


In  Equation  (30)  is  determined  by  the  methods  given  in  Sections  2.14,  2.15 
and  2.19. 

In  the  design  or  project  stage,  it  is  difficult  to  asoribe  values  to 
P"  i.e.  the  density  at  the  throat.  However,  the  density  relationship  through- 
out the  primary  zone  may  be  written:  - 


(31) 


which  is  based  on  a temperature  relationship  assuming  constant  statio  pres- 
sure. G is  a factor  (0  < G < 1)  depending  upon  the  amount  of  heat  release 

p 

between  the  exit  of  the  swirler  and  the  throat.  The  value  of  — = x say, 

P’ 
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oan  usually  be  fixed  with  a reasonable  accuracy,  and  Equation  (31)  reduces 
to 


1 

p ~ G (x  - 1)  + 1 


(32) 


By  substituting  probable  values  for  G,-  p H is  obtained.  In  praotioe 
it  is  doubtful  if  G will  exoeed  0.5  and  generally  0.25  < G<  0.5.  In  the 
previous  analysis  9 talcing  values  of  G of  0.25  and  0.5  varied  the  primary 
loss  factor  some  30  per  cent  and  the  overall  loss  factor  Bane  7 per  cent. 

Thus  the  value  ascribed  to  G is  not  really  oritioal  in  determining  the  over- 
all loss  of  the  ohanber. 

3.0  Primary  stabiliser  losses 

The  fundamental  prinoiple  of  flame  stabilisation  is  to  reduce  the 
local  velooity  and  effeot  a flow  reversal  by  which  fresh  mixture  is  added  to 
the  piloting  region  to  propagate  combustion.  This  is  achieved  by  two  dis- 
tinct forms  of  piloting  system,  viz.  gutter  and  plain  baffle  type  stabilisers. 
The  former  type  are  used  where  high  velooity  conditions  exist  i.e.  ram  jets, 
reheat  eto.,  and  although  considerable  work  is  being  oarried  out  on  gutters 
few  published  notes  are  available . The  plain  baffle  type  stabilisers 
incorporating  a swirler  are  used  in  the  majority  of  aero  engine  and  indus- 
trial type  chambers. 

3.1  Plain  baffles 


These  plain  baffles  are  of  varying  form  although  they  do  preserve 
some  symmetry  in  design.  To  obtain  the  complete  baffle  loss  the  pressure 
loss  of  the  various  free  area  shapes  (holes,  3ooqps,  eto.),  whioh  constitute 
the  baffle  must  be  determined.  When  air  flows  through  these  various  holes 
the  issuing  free  jets  are  conoidal  in  shape  and  hence  give  rise  to  a dis- 
charge coefficient . If  and  Aq  are  the  overall  discharge  ooeff ioient  and 
total  free  area  of  the  baffle  respectively  and  the  various  components  have 

free  areas  A.j , A2,  Aj  and  discharge  coefficients  Ch  , C^_ , 

Cd3  1 * 


then 


°do 


a-i  + a2  Ca2  + a3  od3 

A-j  + A2  + A3  ..... 


s A Cd 
£ A 


(33) 


Thus  any  shape  or  size  of  baffle  can  be  reduced  to  the  simple  case  of  an 
equivalent  hole  in  a flat  plate.  The  neoessaxy  experimental  values  of  dis- 
charge ooeff ioient  are  taken  from  experimental  results  obtained  at  N.G.T.E. 
and  will  be  published10  in  collected  form  shortly.  Briefly,  the  various 
baffles  were  mounted  in  a test  Beotion  and  the  lose  of  total  pressure 
measured  for  a range  of  velocities . Theoretically1 1 , the  pressure  loss  is 
calculable  providing  the  free  area  of  the  baffle,  the  oross-seotional  area 
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to  whioh  the  flow  expands  (in  terms  of  area,  since  this  is  actually  a diffu- 
sion process)  and  the  discharge  ooeffioient  are  known.  The  former  values 
are  obtained  by  actual  measurement  but  the  discharge  ooeffioient  can  only  be 
determined  experimentally. 

3.2  Variables  affecting  the  pressure  loss  of  plain  baffles 

3.21  Effeot  of  velooity 

1 1 

The  theoretical  curves  show  that  the  non-dimensional  loss  factor 
increases  with  Mach  number  which  for  oonstant  statio  temperature  is  propor- 
tional to  the  velooity.  Pressure  loss  tests  on  various  baffle  shapes  have 
shown  that  the  loss  does  in  fact  increase  with  Mach  number  but  at  a reduced 
rate  of  increase  to  that  predicted.  It  was  thought  that  an  increase  in  the 
discharge  ooeffioient  with  Mach  number  might  account  for  the  discrepancy  and 
this  has  now  been  substantiated  by  independent  experiments^.  Variation  of 
(Ca)  with  vena-oantracta  Mach  number  is  shown  in  Figure  10.  For  most  oonbus- 
tion  chambers  the  change  of  C4  is  small,  but  since  the  loss  is  inversely 
proportional  to  its  effeot  will  be  significant. 

3.22  Effect  of  area  ratio 

For  a given  shape  of  hole  the  discharge  ooeffioient  increases  ini- 
tially almost  as  the  square  of  the  area  ratio  as  shown  in  Figure  11  in  whiah 
the  relative  coefficient  is  plotted  against  area  ratio.  This  curve  is  based 
on  values  obtained  from  Reference  1 3 and  by  experiment  and  is  for  sharp 
edged  circular  orifices.  The  equivalent  curve  for  other  shapes  of  orifice 
will  be  slightly  different . 

3.23  Effeot  of  hole  3ize 

The  effeot  of  using  baffles  containing  a similar  total  area  of  holes 
of  different  size  has  not  resulted  in  any  definite  conclusions  being  reached. 
A large  nunber  of  small  holes  would  be  expeoted  to  give  a higher  loss  on 
aooount  of  the  larger  wetted  area  available  for  friotion.  However,  experi- 
mental results  show  the  converse  to  be  true,  i.e.  the  baffle  having  a small 
nunber  of  large  holes  has  a 2 per  oent  higher  pressure  loss.  It  should  bo 
appreciated  that  the  experimental  error  is  of  this  order  and  also  variation 
in  the  diameter  of  the  holes  has  to  be  extremely  small  to  aooount  for  this 
difference. 


3.24  Effeot  of  hole  shape 

The  shape  of  the  hole  for  a given  free  area  does  affect  the  pressure 
loss  by  variation  in  the  discharge  ooeffioient.  Circular  holes  have  the 
lowest  discharge  ooeffioient  for  a given  free  area.  Square  orifioes  have 
slightly  higher  values  of  and  rectangular  and  elliptical  orifioes  with  high 
values  of  major/minor  axis  ratio  higher  values  still.  Typical  minimum  values 
i.e.  corresponding  to  infinite  area  ratio,  are  given  in  Table  I below. 

TABLE  I 


Type 

Ciroular 

Square 

Rectangular 

Elliptical 

Elliptical 

Axis  ratio 

- 

1 

3:1 

2:1 

4:1 

Cd 

0.60 

0.62 

0.63 

0.62 

0.63 

hydraulic 
mean  depth 

0.282/A 

0.25/A 

0.214/A 

0.26/A 

0.204/A 
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The  vena-contraota  ia  formed  by  the  inward  radial  flows  on  the  up- 
stream face  of  the  baffle  acting  on  the  jet  periphery.  For  a hole  in  the 
centre  of  the  baffle  these  contracting  foroes  are  strongest  when  acting 
normal  to  the  jet  surface.  For  a oiroular  hole  the  foroes  act  normal  to  the 
surface  over  the  entire  periphery  and  produce  the  greatest  contraction. 

Thus  oontraotion  coefficient  (Cp)  increases  as  the  hole  shapes  beoome  "less 
oiroular",  i.e.  elliptical  (2:1;,  square,  rectangular,  etc. 

(Cv),  velooity  coefficient,  represents  the  ratio  of  actual  to 
theoretical  velocity  through  the  hole  and  is  due  to  viscosity  and  boundary 
friction.  Henoe  increasing  the  periphery  of  a hole  for  a given  cross- 
sectional  area  results  in  a decrease  of  (Cy) . For  holes  in  thin  plates  Cy 
tends  to  unity  and  as  periphery  variations  (as  shown  in  Table  I where 

hydraulic  mean  depth  = 2£P-£ff. , are  small,  changes  in  (Cv)  are 

periphery 

negligible.  Since  discharge  coefficient  = Cd.  = C0  • Cy  changes  in  C0  will 
be  the  predominating  factor.  Thus  for  maximum  discharge  through  a given 
area  the  hole  shape  should  be  reotangular  with,  for  example,  an  axis  ratio 
of  4:1.  However,  practical  disadvantages  such  as  corner  stress  concentra- 
tions and  manufacturing  difficulties  may  outweigh  the  advantage  of  the  small 
increase  in  (Cd) . 

For  an  annulus  around  a hemispherical  baffle  mean  values  of  0.9  for 
the  discharge  coefficient  were  obtained. 

"Thumbnail"  sooops  have  a discharge  coefficient  olosely  approaching 

unity. 


3.25  Effect  of  hole  arrangement 

No  general  conclusions  may  be  drawn  from  the  disposition  of  holes  in 
a baffle.  Various  arrangements  of  holes,  for  a constant  area  ratio,  lead  to 
negligible  changes  in  the  overall  loss  faotor. 

3.26  Effect  of  hole  inclination 


To  determine  the  effect  of  inclination  of  the  plane  of  the  hole  to 
the  air  stream  a series  of  cones  were  tested  in  vMoh  the  cone  angle  was 
varied  but  the  area  ratio  and  hole  arrangement  remained  the  same.  When  the 
holes  were  placed  normal  to  the  airstream  minimum  loss  was  obtained.  As  the 
angle  between  the  axes  of  the  holes  and  the  airstream  (6)  increased  the  loss 
increased  approximately  as  cos^  6 as  shewn  in  Figure  12.  This  is  to  be 
expeoted  sinoe  the  projected  area  of  the  holes  on  a plane  normal  to  the  air 
flow  is  direotly  proportional  to  oos  6 and  loss  is  proportional  to  the  square 
of  the  area  ratio.  Figure  13  shows  relative  loss  defined  as 

Jag£  iAP^.ation  6 plotted  against  6.  Placing  the  cone 

loss  faotor  at  0 = 0 (i.e.  minimum  loss; 

apex  upstream  or  downstream  had  no  measurable  effeot  on  the  loss. 


3.27  Effeot  of  turbulenoe 

Reference  14  gives  details  of  experiments  carried  out  cm  a series  of 
flat  plates  whioh  illustrate  the  effeot  of  turbulenoe  on  drag.  Figure  14 
shows  the  variation  of  pressure  drop  coefficient  (statio  pressure  difference/ 
free  stream  velooity  head)  with  percentage  turbulenoe.  The  percentage  turbu- 
lence is  defined  as 


root  mean  square  of  speed  fluctuation  ^qq 
average  speed 
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The  turbulenoe  level  was  varied  by  placing  large  wire  diameter,  large 
mesh  gauzes  upstream  of  the  test  section.  Considering  praotioal  applications, 
the  change  of  percentage  turbulenoe  is  usually  small  in  a given  test  set  up, 
but  this  feature  of  drag  increase  with  percentage  turbulenoe  is  important 
when  comparing  pressure  loss  measurements  made  on  an  identioal  component  on 
two  dissimilar  rigs.  However,  reference  to  Figure  14  shows  that  percentage 
turbulenoe  changes  will  only  account  for  small  differences  in  pressure  loss. 

3.3  Gutter  stabilisers 


The  loss  due  to  gutters  is  mainly  an  expansion  loss  arising  fran  the 
fuel  injeotor  situated  in  the  high  velooity  throat  and  also  the  diffusion 
loss  up  to  the  ohanber  oross-section  fran  the  downstream  end  of  the  gutter. 
For  incompressible  flow  the  loss  is  (X  - 1)2  and  includes  a discharge 
ooeffioient  far  the  gutter.  For  inoluded  gutter  angles  up  to  15°  the  value 
of  C<3.  is  about  unity.  For  higher  angles  the  Cd  decreases  fairly  rapidly, 
probably  following  a oosine  law,  but  this  is  merely  a hypothesis  whioh,  - 
although  qualitatively  correct,  should  be  oonfirmed  experimentally  before 
being  used  indiscriminat e iy . If  the  throat  velocity  is  greater  than 
200  ft./seo.  the  curves  of  Reference  11  should  be  used  to  allow  for  compres- 
sibility in  determining  the  pressure  loss. 


For  hot  running  the  fundamental  pressure  loss  due  to  heat  addition 
(see  Seotion  6.0)  is  added  to  the  oold  loss.  The  result  obtained  may  be  high 
compared  with  the  experimental  value.  This  is  due  to  the  aerodynamic  flow 
pattern  around  the  gutter  being  signifioantly  altered  by  combustion.  The 
principal  effects  of  combustion  are  to  reduce  the  strength  of  the  reverse 
flow  (and  henoe  the  pressure  loss)  and  to  increase  the  length  and  breadth  of 
the  wake.  A further  contribution  to  the  loss  factor  is  the  dissipation  of 
the  upstream  component  of  the  fuel  momentum  when  injected  in  the  throat.  If 
the  inlet  air  and  fuel  tenperatures  are  substantially  the  same,  inareas e in 
fuel  flow  results  in  an  increase  in  pressure  loss  (of  the  order  3**5  per  oent), 
but  if  the  air  temperature  is  high  compared  with  the  fuel  the  pressure  loss 
tends  to  decrease.  This  later  phenomenon  is  due  to  the  reduction  in  air 
temperature  due  to  fuel  vaporisation.  The  presenoe  of  the  fuel  increases  the 
effective  blockage  at  the  throat,  and  3inoe  the  throat  velooity  and  permanent 
blookage  are  both  high,  exerts  a measurable  effeot  on  the  loss.  If  the 
throat  section  is  long  friction  effeots  must  be  taken  into  account  by  the 
modified  "Fanning  Equation" 


iiM  = a £ 

dt  . D 


(34) 


for  reotangular  or  annular  oross-sectional  areas  the  equivalent  diameter 
(de)  is  used,  (f)  will  vary  between  0.002  and  0.008  depending  on  the 
Reynolds  number  as  shown  in  Figure  15. 

The  effeot  on  pressure  loss  of  using  skirted  gutters  (see  Figure  16, 
(a)  and  (b))  as  opposed  to  the  conventional  type  is  negligible,  although  an 
improvement  in  flame  stability  may  result.  The  use  of  "finger"  type  flame 
spreaders  attached  to  the  downstream  end  of  the  gutter  gives  rise  to  a snail 
increase  in  the  loss  whioh  is  accounted  for  approximately  by  the  loss  due  to 
flow  through  the  projected  free  area  of  the  fingers  in  the  plane  of  the 
gutter  base  as  shown  in  Figure  l6(o) . This  loss  will  probably  be  a little 
higher  than  the  more  gradual  loss  occurring  along  the  fingers,  but  does  give 
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a 'basis  f or  analytical  determination.  Diffusion  losses  oan  be  treated  by 
the  method  given  in  Section  7.1  • 

From  the  preceding  paragraph  it  is  obvious  that  the  pressure  loss 
pioture  is  far  from  complete,  but  correlation  of  the  results  of  many  experi- 
ments now  in  progress  will  improve  the  position. 

4,0  Pooling  losses 

The  main  types  of  oooling  device  in  use  at  the  present  time  are  the 
louvre,  porous  wall  and  boundary  layer  systems.  For  a detailed  analysis  and 
description  of  these  systems  Reference  1 5 should  be  consulted.  From  the 
point  of  view  of  pressure  loss  no  new  problems  are  involved,  each  system 
merely  utilising  the  available  pressure  difference  between  the  primary  and 
secondary  flew  paths. 

4.1  PorouB  wall 


This  method  of  oooling  is  among  the  more  efficient  and  is  amenable  to 
analytical  treatment.  The  pressure  drop  for  laminar  flow  through  a porous 
material  is  given  by  D’Aroy’s  equation 


Q.L.  _ _Z_  p2 
Pi2 *  - P22  144  2 p2  H 


(35) 


2 

where  p^  and  pg  are  the  air  pressures  in  lb, /ft.  on  either  Bide  of  the 
porous  wall  and  Z is  the  coefficient  of  permeability  and  has  dimensions  of 
an  area,  usually  square  inches.  For  the  small  pressure  differences  available 
in  ccrrbustion  chambers 


2 2 

Pi  - Pg  n 2 P2  • Ap  where  Ap  = pressure  drop 


thus 


Q.L.  Z p 
Ap  ~ 144  M 


(36) 


However,  Z must  be  determined  experimentally  in  the  first  instance, 
and  may  decrease  with  operating  time  due  to  deposition  in  the  pores. 

Typical  values  of  Z are  1 0“® to  1 0"1 0 in2, 

4,2  "Louvred"  surface  oooling 

The  "louvred"  wall  is  essentially  a mode  cf  oonstruoticn  (British 
Patent  No. 642, 257  held  by  "Shell"  Refining  and  Marketing  Company  Limited) 
by  whioh  the  effective  area  for  heat  transfer  is  considerably  increased. 

The  surface  to  be  oooled  is  constructed  so  that  there  are  many  small 
independent  passages  along  whioh  the  oooling  air  may  flow  radially,  finally 
emerging  to  mix  with  the  primary  stream.  To  estimate  the  pressure  drop 
associated  with  the  flow  of  air  through  the  passages  in  the  "louvred  wall" 
under  turbulent  flow  conditions,  Blaxiua'  Equation  is  usedt- 
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„ 2 

- Qlj16  . Pv  L 

Re0 *25  2 de 


(37) 


for  laminar  flew  conditions: - 


AP  = 


96 

Re 


2 

£V_  _L_ 

2 * de 


(38) 


The  oriterion  for  turbulent  or  laminar  flow  is  whether  Re  is  above  or 

1 2 

below  2,000.  In  addition  the  injection  loss  a - p V should  be  added  to 
either  Equation  (37)  or  (38)  to  give  the  oomplete  loss. 


4-3  Combination  of  external,  air  flow  and  localised 
air  injection  pooling 

This  method  requires  the  cooling  air  to  flow  in  an  annular  sheath  in 
an  upstream  axial  direction  and  then  to  inject  it  through  small  holes  in  the 
flame  tube  into  the  high  temperature  side  where  it  forms  a blanketing  annular 
layer.  The  pressure  drop  i3  the  sum  of  the  friction  drop  given  by  Equations 
(37)  or  (38)  and  the  injection  loss  whioh  will  be  due  to  accelerating  the 
air  up  to  the  required  injection  velocity.  The  latter  loss  is  given  by 


where  is  the  velocity  based  on  total  port  area  and  0.61  is  the  discharge 
coefficient . 

The  overall  loss  for  this  cooling  system  is 

fl?.  L , P / Vh  \ 

Re0*25  2 de  2 '0.61' 


5.0  Mixing  losses 

Up  to  this  Seotion  most  of  the  information  is  complete  and  valid  for 
all  types  of  combustion  chamber  but  an  incomplete  knowledge  of  the  mixing 
process  restricts  the  application  to  low  speed  chanters. 

The  pressure  loss  due  to  mixing  is  probably  the  most  difficult  loss 
to  determine  analytically  without  some  experimental  assistance,  since  it 
affeots  both  the  primary  and  secondary  streams.  The  part  of  the  mixing  loss 
attributable  to  the  seoohdary  oirouit  is  almost  entirely  due  to  expansion 
through  the  mixing  holes.  The  loss  associated  with  the  primary  oirouit  is 
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made  up  of  the  flow  through  the  effective  "blockage  due  to  the  radial  "Bpokes" 
of  oold  air  and  the  subsequent  m&oro-turbulenoe. 

5.1  Secondary  mixing  losses 

As  stated  in  the  previous  Seotion  the  seocndaiy  mixing  loss  is  given 
approximately' by  the  velocity  head  through  the  holes.  This  requires  a 
knowledge  of  the  discharge  coefficient,  which  is  subject  to  a wide  variation 
depending  on  hole  area,  outer  duot  area  and  the  percentage  flow  from  the 

outer  duct  through  the  hole.  Figure  17  is  a curve  of  versus  a factor  2 

where  F is  the  percentage  flow  from  the  outer  duot  and  B is  the  ratio  of 
hole  area/outer  duct  area.  This  curve  was  taken  from  Reference  1 6 and  is 
the  result  of  water  model  tests  with  hole  sizes  ranging  from  0.6  to  1 .9  in. 
diameter.  It  is  satisfactory  to  determine  the  percentage  flow  through  the 
hole  on  an  area  basis. 

1 7 

Darling  ' has  also  studied  this  problem  using  air  as  the  flow  medium 
and  presents  his  values  of  discharge  ooeffioient  as  a function  of  the 

"Approach  Velooity  Factor",  i.e.  ~ where  V-j  is  the  mean  velocity  in  the 

approach  channel,  and  V2  is  the  mean  velooity  through  the  hole.  The  nunfcer 
of  experimental  points  taken  are  less  than  in  Reference  17  and  only  one 
size  of  hole  was  used.  Darling's  results  have  been  plotted  on  the  same 
absoissa  as  the  Lucas  results  in  Figure  17.  The  curves  are  of  similar  shape 
although  the  curve  for  air  is  some  7 per  cent  higher.  For  equal  conditions 
of  flow  the  discharge  ooeffioient  for  air  would  be  higher  due  to  compressi- 
bility although  by  a very  small  amount.  The  real  difference  appears  to  be 
due  to  the  positioning  of  the  static  taps  on  the  two  separate  rigs.  For  the 
water  model  they  are  situated  in  the  annulus  some  2f  in.  upstream  of  the 
injection  hole  axis  whereas  for  the  air  tests  the  tap  was  situated  on  the 
outer  annulus  wall  directly  above  the  centre  of  the  hole.  The  maximum  value 
of  C4  obtained  in  Referenoe  17  is  higher  than  anticipated  for  this  type  of 
discharge . The  true  values  for  air  are  probably  a little  higher  than  the 
water  results  although  negligible  error  will  result  in  applying  these 
directly  to  air  calculations.* 

The  secondary  pressure  loss  due  to  mixing  will  then  be  given  by 


Ca  being  obtained  from  Figure  1 7 . 


*This  statement  is  oonfirmed  by  an  Amerioan  Report  "Can  Burner  Hole  Discharge 
Ooeffioient  Investigation"  Consolidated  Vultee  Aircraft  Corporation  No. 6149, 
just  reoeived. 
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Far  holes  inclined  to  the  direction  of  fleer  the  discharge  coefficient 
obtained  from  Figure  17  should  be  increased  by  the  root  of  the  relative  loss 

factor  since  C^oe^i  . For  example,  if  the  mixer  has  a semi-angle  of  15° 

then  with  reference  to  Figure  13,  0 a 75°,  and  Ca  obtained  from  Figure  16  is 

multiplied  by  lljh§£  . 

%/ 1.452 

The  preceding  statements  assume  that  the  hot  stream  effects  are 
negligible.  This  is  probably  true  for  low  speed  industrial  type  chambers 
but  evidence  from  experiments  now  in  progress  suggests  that  the  hot  stream 
momentum  substantially  affects  the  result  and  reduces  the  value  of  the  pres- 
sure drop  as  given  by  Equation  (41 ) . 

5.2  Primary  mixing  losses 

Losses  in  the  hot  stream  from  the  inject  ion  plane  to  the  "mixed”  plane 
are  approximately  half  the  velocity  head  at  the  plane  of  injection  and  are 
thus  very  small.  For  very  large  or  industrial  type  ahambera  it  oan  be 
regarded  as  negligible.  This  part  of  the  work  will  be  in  a much  mare  exact 
form  when  the  results  of  mixing  experiments  new  in  progress  are  available. 

6.0  Heat  addition  losses 


If,  as  is  usual,  the  oonbustion  occurs  in  a parallel  duct  immediately 
downstream  of  the  primary  baffle  the  "fundamental"  loss  of  pressure  is  given 
by 


AP.ip,  v,2(|l-l)  <«> 

and  if  the  static  pressure  difference  is  small 

1 2 / To  \ 

AP  »-  Pi  V1  ^ - 1 J (42a) 


In  the  oase  of  a varying  oross-seotional  area  in  the  flame  tube,  it  is 
best  to  consider  in  detail  the  relative  proportions  of  heat  release  as  in 
Section  2.3. 

7.0  Miscellaneous  losses 

7.1  Diffusion  losses 

For  various  reasons  the  reduction  of  velooity  in  the  compressor 
diffuser  is  often  limited  and  the  inlet  velocity  to  the  oonbustion  chamber  is 
frequently  greater  than  300  ft  ./sec.,  the  exact  value  depending  to  a large 
extent  on  the  type  of  oenpreasor.  Typical  values  for  the  velooity  in  the 
seoendazy  annulus  are  of  the  order  of  150  ft  ,/seo.  and  it  is  neoesssry  to 
reduoe  the  inlet  air  velooity  to  that  existing  in  the  aeoondary  annulus  as 
efficiently  as  possible.  The  efficiency  of  a diffuser  may  be  defined  by  a 
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factor  (e)  which  givas  the  effioienoy  of  conversion  of  velocity  head,  to 
static  pressure 


»2-*1 


(43) 


. . total  pressure  loss  P.,  - Pg  « (1 


e)  -L  (v^  - v22)  ..  ..  (44) 


loss  factor  4 = (1 


(45) 


Values  of  ( e)  have  been  taken  from  Reference  1 8 which  agree  with 
experimental  results  given  in  Reference  19  and  are  plotted  in  Figure  18 
against  total  diffuser  angle  8.  A recent  report20  has  shown  that  assyro- 
metry  of  the  inlet  velocity  distribution  has  a marked  effect  on  diffuser 
effioienoy  especially  for  large  diffuser  angles.  A low  velooity  region  near 
the  wall  is  equally  undesirable. 


7.2  Losses  due  to  bends 


Although  not  explicitly  a component  .of  the  combustion  chamber,  bend 
entries  and  exits  for  oorbustion  chambers  are  relatively  common  and  their 
loss  is  frequently  inoluded  in  the  overall  ahamber  loss  figure.  Aoourate 
data  far  the  losses  in  bends  is  given  in  Reference  21 , but  in  general  terns 
it  can  be  stated  that,  for  a bend  without  diffusion  and  with  a directional 
ohange  not  exceeding  90°,  and  having  a mean  radius  not  less  than  1 .5  times 
the  duct  diameter  or  passage  width,  the  pressure  loss  will  not  exceed  half 
the  velooity  head.  The  loss  round  a sharp  bend  can  be  reduced  by  imparting 
an  acceleration  to  the  air. 

Cascade  bends  are  now  universally  employed  in  gas  turbine  systems  by 
virtue  of  their  effioienoy  both  in  terms  of  pressure  drop  and  their  ability 
to  turn  the  air  through  a desired  angle.  Reference  22  gives  the  design 
details  and  procedure  for  constructing  a bend  in  which  the  blades  are  spaced 
in  an  arithmetic  progression  from  the  inside  radius.  The  pressure  loss 
associated  with  such  a bend  is  affected  by  size  and  manufacturing  variations 
(especially  internal  finish)  but  a loss  figure  of  25  per  oent  of  the  velooity 
head  through  the  bend  is  sufficiently  accurate  for  most  purposes. 

7.3  Losses  due  to  oorrugated  spacers 

This  form  of  construction  is  now  used  frequently  as  a mechanical 
spacer  for  skin  pooling  of  oonfcustion  chamber  walls.  The  discharge  coeffi- 
cient cf  this  spacer  was  investigated?*  on  a water  model  and  found  to  be 
0.8  when  based  on  the  drawing  dimensions  and  0.9  in  terns  of  the  actual 
measured  areas.  The  variation  in  drawing  and  measured  dimensions  is  due  to 
manufacturing  difficulties  principally  in  the  welding  operation.  For  design 
purposes  the  estimated  area  of  the  section  is  used  for  which  Cd  equals  0.8. 
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7 A Friotlon  loaaea 


In  the  majority  of  ohambers  the  friction  losses  are  usually  negligible 
compared  with  the  individual  ocmponent  losses  but  in  a few  isolated  oases 
there  are  long  seotions  where  friotional  affeots  are  measurable . 

The  pressure  drop  is  given  by  the  modified  "Fanning  Equation" 


iM 

04 


(3V) 


f being  obtained  from  Figure  1 5 . 

For  irregular  shaped  ducts  and  annuli  the  hydraulio  mean  diameter  d@ 
is  used  for  D. 


A 

i.e.  D = 4 g 


8.0  Overall  chamber  loss 

Having  considered  in  acme  detail  the  pressure  losses  caused  by  the 
various  components  of  the  chamber,  it  is  now  necessary  to  see  how  they  may 
be  linked  to  give  a value  of  the  loss  coefficient  for  a particular  flow  path. 
For  consecutive  losses  in  a flow  path  the  overall  loss  coefficient  is  merely 
the  arithmetic  sum  of  the  individual  loss  factors  provided  they  are  expressed 
in  terms  of  the  same  reference  velocity  head.  For  losses  occurring  in  paral- 
lel oircuits  the  method  of  Probert  and  Kielland1  iB  used.  A loss  coefficient 
is  applied  to  each  flow  path  such  that  the  total  head  loss  in  the  stream  is 
equal  to  the  loss  coefficient  times  the  velooity  head  at  some  reference  area. 
On  the  further  assumption  that  the  static  pressures  are  equal  in  both  streams 
at  divergence  and  confluence  an  expression  for  the  overall  loss  factor  is 
obtained.  While  this  method  proves  satisfactory  for  the  simpler  types  of 
chamber  a considerable  amount  of  calculation  is  required  if  there  are  more 
than  two  general  flow  paths.  Also,  because  of  the  "step-by-step"  method  of 
calculation,  if  the  loss  factor  of  one  of  the  components  is  changed  a complete 
recalculation  is  neoessary. 

Reference  2 is  based  on  the  same  principles  and  assumptions  as  stated 
above  but  as  shown  in  Appendix  II  reduces  the  complexity  and  quantity  of 
calculation. 


If  fa  s pressure  loss  faotor  of  a 

oirouit  in  terms  of  velooity  • 
head  at  area  x 

and  fa  a same  loss  in  terms  of  velocity 
head  at  area  y 

then  ••  ••  ..  (46) 

fa  \y/ 
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and  in  the  event  of  a density  ohange 


f = (f ) (If) 


Thus  by  applying  this  relation  it  is  possible  to  express  all  indivi- 
dual loss  factors  in  terns  of  the  velooity  heads,  due  to  each  flow,  at  the 
same  oross-seotional  area.  This  referenoe  area  is  purely  arbitrary  and  oan 
be  the  ahamber  entry  area,  flame  tube  area  or  oaaing  area. 

Since  the  total  head  drop  of  any  flew  circuit  in  the  chamber  must  be 
the  same 


<1|  a ^2  = ^3  = * •tl*  **  ••  *•  ••  (48) 


where  q.j , q2  eta . are  the  velooity  heads  due  to  the  individual  flows 

in  the  referenoe  area,  and  etc.  are  the  loss  factors  expressed 

in  terms  of  the  velooity  head  at  the  reference  area  by  means  of  Equation 

(47). 

2 

But  qa  W since  p is  assumed  constant  at  the  reference  area  for  all 

flows 

Wi  = W 2Jfa  = etc (49) 

and  the  overall  loss  factor  by 

* = A (“)  = h (—)  eto (50) 


also  sinoe  the  sum  of  the  percentage  flows  through  each  circuit  must  equal 
the  total  flow 


W = 100  « w1  + W2  + IflTj  eto (51) 

thus  any  required  circuit  flow  say  W-j  is  given  by: 


w1  = 100  - w2  - 

= 100  - w>  (£)*  - *'  (rf 


assuming  there  is  a total  of  three  circuits. 
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9.0  Conclusions 


By  means  of  the  analysis  of  component  pressure  losses  in  this  Report 
it  should  be  possible  to  make  a reasonably  accurate  theoretical  calculation 
of  the  oold  air  flow  distribution  and  overall  loss  factor  of  a combustion 
chamber.  Certain  limitations  in  our  knowledge  of  compressible  flow  charac- 
teristics especially  mixing  of  gas  streams,  imposes  a restriction  an  the 
aoouraoy  for  high  velocity  chanters . This  contingency  will  be  obviated  by 
experimental  work  now  in  hand.  The  comparison  between  calculated  and 
measured  pressure  drop  for  a typical  combustion  chanter  as  shown  in  Append 
dix  II  is  good.  The  percentage  difference  may  be  fortuitous  but  the  prospeots 
of  calculating  the  oold  pressure  drop  of  a ohamber  from  the  design  drawing 
with  an  aoouraoy  of  + 5 per  cent  seems  favourable.  Assuming  the  mixing 
experiments  improve  the  'hot'  pressure  loss  calculations,  the  method  can 
probably  be  further  refined  by  comparing  oaloulated  and  measured  results 
from  a variety  of  chambers. 
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APPENDIX  I 
S.ynfrols 

Cross-seotional  area 

Ratio hole  ■££ — (see  Figure  17) 

outer  duot  area 

Constant  for  vortex  decay  law 

Blade  chord 

Discharge  coefficient 

Lift  coefficient  (see  Equations  9 and  10) 

Outer  diameter 

Inner  diameter 

Equivalent  diameter  = 4 x 

Diffuser  efficiency 


- ft. 

dimensionless 


cross-sectional  area 
perimeter 


Percentage  flow  frc»n  cuter  duot  (see  Figure  17) 
Friotian  factor  (see  Equation  34) 

Heat  release  factor  (see  Equation  31) 

Total  energy  per  unit  mass 
Secondary  loss  factor  (3ee  Equation  10) 

Length 
Mach  number 

Mach  number  at  vena-contraota 
d2 

Area  ratio  = — s 
v 

Index  in  vortex  decay  law 

Total  pressure 

Statio  pressure 

Total  pressure  loss 

Static  pressure  loss 

Mass  flow  per  unit  ooal ed  surface  area 


dimensionless 

dimensionless 

- ft. 

- ft. 

- ft. 

dimensionless 

dimensionless 

dimensionless 

dimensionless 

- ft.  sec.  * 
dimensionless 

- ft. 

dimensionless 

dimensionless 

dimensionless 

dimensionless 

-2 

- Ib.  ft. 

- lb.  ft."2 


lb.  ft. 


-2 


-2 


- lb.  ft. 


- slugs. 

seoT”*  ft.”2 


Radii 


- ft. 
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Reynolds  number 

Perimeter 

Pitch 

Blade  thiokness 
Absolute  velocity 
Velocity  through  hole 
Whirl  velooity 
Weight  flow 
Area 

Profile  loss  coefficient 
Secondary  loss  coefficient 
Total  loss  coefficient 
Area 

Coefficient  of  permeability  (see  Equation  36) 

Outlet  air  angle 

Blade  outlet  angle 

Ratio  of  specific  heats 

Baffle  semi-cone  angle 

Effective  area  ratio 

Viaoosity 

Density 

Flame  tube  density  (see  Figure  9) 

Throat  density  (see  Figure  9) 


dimensionless 

- ft. 

- ft. 

- ft. 

- ft.  sec.”1 

- ft.  sec.”1 

-1 

- ft.  860. 

-1 

- lb.  £30. 


dimensionless 

dimensionless 

dimensionless 


2 

~ m. 


dimen iionloss 
dimensionless 

dimensionless 

dimensionless 

dimensionless 

- slugs. 
ft.-1  sec.’1 

- slugs .ft 

- slugs.ft.-^ 

- slugs. ft.“^ 


Pitch  ohord  ratio  = £ 


Loss  ooeffioient 


Pi  - ?2 


dimonsi  .enless 
dimensionless 


Swirl  angle  for  vortex  swirl or  (see  Section  2. 26) dimensionless 
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Suffioea 


( )o 

a 

Known  condition  usually  inner  radius 

( )1 

a 

Entry  or  initial  condition 

<>2 

a 

Outlet  or  final  condition 

( )p 

a 

Pertaining  to  flame  tube 

<>. 

» 

Pertaining  to  mean  radius 

( )s 

a 

Pertaining  to  swirler  or  secondary 

<>„ 

= 

Throat  condition 
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Air  flow  distribution  and  overall  loss  f apt or 


for  a conventional  ohagfaer 


(Rolls  Royoe  R.M.60  Model) 


Aa  oan  be  seen  from  Figure  the  air  flow  is  divided  into  eigiht 
separate  flow  circuits . Eaoh  individual  loss  factor  will  be  expressed  in 
terns  of  the  velocity  head  pertaining  to  the  overall  ohamber  cross-sectional 
area. 


To  determine  the  "hot"  distribution  case  at  a given  temperature  ratio 
the  oold  distribution  is  used  to  oaloulate  the  primary  oombustion  zone 
temperature.  Striotly,  a method  of  successive  approximation  should  be  used 
to  allow  for  small  redistributions  of  air  flow  but  the  magnitude  of  the 
errors  involved  and  the  general  accuracy  of  the  method  as  a whole  do  not 
warrant  it . 

Calculation  of  individual  loss  factors 

(1 ) Expansion  ratio  through  primary  orifice 


= !a2L  = 3.47  i.e.  m = 0.312 
0.84 


From  Figure  1 1 


Cd  = 0.6  x 1 .058  = 0.635 


The  effect  of  the  shoulder  will  certainly  reduce  the  discharge  and  a 
Ca.  °f  0.6  is  used. 


2 

Loss  through  orifice  = (A  - 1 ) - 


- 22.8 


I 

In  terms  of  reference  /i 


22.8 


i±222 


Considering  the  swirler 


54 


Ap  a 21 .3  in. 


Ag  a 2.6  in. 


By  Equation  (12) 


*F 


1 

(0.5878)2 


1 


8 222 
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In  terms  of  reference  area  / = 222  x 1 ■ 


= 725 

The  overall  loss  for  the  two  resistances  in  series  is  the  algebraio  sum  of 
the  loss  factors  (when  expressed  in  terms  of  the  same  area) . 


. . Loss  through  No.1  Circuit  -fa  =4,715 
(2)  Loss  through  corrugated  spacer 

2 

Free  area  = 1 .07  in. 

2 

Expanded  area  = 2.43  in. 
from  Section  7.30  C<j  = 0.8 


.2  .2 

Loss  in  terms  of  reference  area  = [——ihL—-  - 1 j ( ) = 850 

\1 .07  x 0.8  / \2.43/ 


'Expansion'  I033  after  spacer  in  terms  of  reference  are 

2 . .2 


■{(lid)  - '}  (M)-’” 


Loss  through  No. 2 Circuit  = fa  - 1 ,047 
(3)  Loss  through  primary  holes 

Firstly,  the  C<j  of  the  holes  must  be  estimated  by  the  method  outlined 
in  Section  5.10  and  Figure  17. 

F is  determined  on  an  area  basis  only 


100  x 1.39 

1 .39  + 0.48  + 0.55  + 3.49  + 1.99  + 3.49 


139 

11.39 


12,2  per  cent 
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B = ij-22  = 0.0837 

16.6 

P - 12*2  = 146 


B 0.0837 
From  Figure  17  C^  = 0.582 

2 


Loss  in  terms 


of  reference  area  = ( — 1 — _ -A  ( — 

\1  .39  x 0.582  } \21 .3/ 

^3  = 2,090 

(4)  Loss  through  first  row  of  cooling  holes 


w - 

100  x 0.48 

10.0 

B = 

0.48 

16.6 

F _ 

4.80 

B = 

0.0289 

4.80  per  oent 
0.0289 
166 


From  Figure  17  Cd  = 0.595 


Sinoe  these  holes  are  inclined  at  an  angle  of  17°  the  discharge 
ooeffioient  is  increased  (see  Section  5.10  and  Figure  13). 


Cd  = 0.595  x l^2  = 0.603 
J 1.442 


It  is  assumed  that  the  air  entering  these  holes  forms  an  annular  sheath 
which  does  not  substantially  increase  in  thickness  as  it  flows  downstream. 


loss  factor  /Sl  =»  ( 

\ 0.603  x 0.48 


^ = 13,800 
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(6)  Losa  through  first  row  of  mixing  holes 

F = 19P-i ll/it?.  = 38.9  per  oent 
8.975 

B = 1^2  = 0.309 

11.3 

I.  &SL  =126 

B 0.309 

. . Prom  Figure  17  C4  a 0.56k 
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Mixing  loss  «*  — P ^ 

2 VV 


fa  ■ 383 


Loss  through  third  set  of  cooling  holes 
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F = 


100  x 1 .99  _ 1 99 
1.99  + 3.49  5.48 


= 36.3  per  cent 


b = hSl 
11 .3 


0.176 


F _ 36.3 

B “ 0.176 


206 


1 

) 


. . Cd  0.608 

Since  inclination  is  20°  Cd  is  further  increased.  From  Figure  13 


Cd  = 0.608  x = 0.62 

V .426 

By  Equation  (41 ) the  loss  factor  in  terms  of  the  hole 

*rea-  (0^2)  ■ 2-6 

<4  - 2-6  (i?i)2  - 971 

fa  = 974 


I 
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Equation  (41 ) was  used  as  it  is  very  difficult  to  decide  to  which  eff eotive 
area  the  injected  air  eventually  "expands". 

(8)  Loss  through  final  mixing  holes 

P = 100  per  oent 


Cold  air  distribution 

h. 

H 

K 

H 

h 

h 

h 


£ = 1QQ  = 238 
B 0.42 


Cd 

= 0.61 

2 

2 

y8  = (o7 

\ / 38.5' 
/ \ 3.49< 

) = 327 

r-  f 

CM 

II 

XP 

= 4,715 

II 

T— 

"Q. 

*> 

68.7 

= 1,047 

V ^2  ® 

32.4 

= 2,090 

45.7 

= 13,800 

11 

* 

•> 

117.6 

= 10,700 

V ^5  = 

107.0 

W\ 

$ 

II 

= 

19.6 

= 974 

^7  = 

31 .2 

- 327 

II 

18.1 

I 
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By  Equation  (51) 


100 


+ 68^1  + 68j_  

32.1,.  45.7  ' 117.6  ' 107.0  ' 19.6  ' 31.2  ' 18.1 


, 68.7  , 68.7  , 68. 


100 

1 + 2.12  + 1.50  + 0.58  + 0.64  + 3.51  + 2.20  + 3.80 


100 

15.35 


^ a 6.5  per  oent 


F 2 


100 


100 100 

0.47  + 1 + 0.71  + 0.27  + 0.30  + 1 .65  + 1 .04  + 1 .79  ° 7.23 


W2  = 13.8  per  oent 


W,  = 


100 


3 4t5uZ  + + i + Jt^2  + + + + 

68.7  32.4  117.6  1 07.0  1 9.6  31.2  1 8.1 


100 100 

0.66  + 1 .42  + 1 + 0.39  + 0.43  + 2.34  + 1 .47  + 2.53  10.24 


Wj  a 9*8  per  oent 


100 

4 = HLi  + Uii  + nz^  + i + HM  + ilLi  + ULi  + 

68.7  32.4  45.7  107  19.6  31.2  18.1 


100 100 

1 .71  ♦ 3.63  ♦ 2.57  + 1 + 1 .10  + 6.00  + 3.77  + 6.50  “ 26.28 

3.8  per  oent 
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12Z_  + 12L  + 12L.  + _12Z_  + i + 12L.  + i£L  + 121. 

68.7  32.4  45.7  117.6  19.6  31.2  18.1 

100 100 

1.56  + 3.30  + 2.34  + 0.91  + 1 + 5.42  + 3.43  + 5.90  " 23.86 


W5  » 4.2  per  oent 


* w 

"6  = 1M  * 13A  + 12.6  + 7ISJ777I 2*  V12A 

68.7  32.4  45.7  117.6  107  31.2  18.1 

100 100 

0.28  + 0.60  + 0.43  + 0.17  + 0.18  + 1 + O.63  + 1 .08  " 4.37 


Wg  = 22.8  per  oent 


7 = 2L4  + 21A  + il^2  + -2^  + _J1*2  + + 1 +iuz 

68.7  32.4  45.7  117.6  107.0  19.6  18.1 

100 100 

" 0.45  + 0.96  + 0.68  + 0.26  + O.29  + 1.59  + 1 + 1 .72  “ 6.95 


W7  = 14.4  per  oent 


"8  = 1M  + IM  + 18*1  + 18.1  + 18.1  + 18*1  + Jlii  + 1 

68.7  32.4  45.7  117.6  1 07.0  1 9.6  31.2 

100 100 

= 0.26  + 0.56  + 0.40  + 0.15  + 0.17  + 0.92  + 0.58  + 1 3 4.04 


Wg  = 24.7  per  oent 


Cheok:-  6.5  + 13.8  + 9.8  + 3.8  ♦ 4.2  + 22.8  + 14.4  + 24.7  » 100  per  oent 
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i Overall  oold  loss  factor  by  Equation  (50) 

t 


= 19.9  in  terms  of  referenoe  velocity  beads 


2 2 
. /W2\  ,13. 8\ 


The  measured  value  of  the  oold  pressure  loss  factor  was  20.7  an  error 
of  about  4-  par  oent. 


Hot  pressure  loss 

To  determine  the  effect  of  heat  addition  it  is  neoessary  to  arrive  at 
a value  for  the  primary  temperature . 


Using  the  previously  determined  air  flow  distribution  and  assuming 
oirouits  1,2  and  3 constitute  the  primary  air  flow. 


Percentage  primary  air  =»  6.5  + 13.8  + 9.8  * 30.1  per  oent 

Neglecting  speoifio  heat  variation  and  assuming:  - 

Inlet  temperature  = 200°C. 

Outlet  temperature  = 700°C. 

If  is  the  primary  absolute  temperature 

then  30.1  T-i  + 69.9.4.73  = 100.973 

T 97.300  - 33,000  64,300 

* 1 " 30.1  “ 30.1 

= 2,130°K. 

By  Equation  (42a) 

Heat  addition  loss  factor  ■ ( - 1 \ 

\ 473  / 
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and  in  terms  of  referenoe  area 


= -HA 

V 1.73  / \21 .3/ 


Overall  primary  loss  factor  excluding  oombustion  loss  is  by  Equation 


. / W \ ,100  \ 

K = *( ) = 19.9  X ( ) = 

P \W-|  + W2  + W3/  \30.1  / 


a 220 


New  primary  los3  factor  including  oombustion  loss  will  be  220  + 11  A 
231  A. 

Assuming  the  secondary  loss  factor  remains  constant 


(0.699)' 


. the  new  distribution  is 


SJ231'4  = Qs>0-8 


peroentage  through  primary 


= 29.5  per  cent 


231  A 


Thus  heat  addition  has  reduced  the  primary  total  flc«r  by  30.1  - 29.5 
= 0.6  per  pent. 


The  new  hot  loss  faotor  = 231  A (0.295)  = 20.2 


The  measured  hot  loss  faotor  for  the  assumed  temperature  rise  was 
25,  an  error  of  about  20  per  cent. 


CCNPIDBNTIAIi 


The  equation  for  radial  equilibrium  in  vortex  flow  is 

P dr  r 

The  total  energy  unit/mass  at  any  radius  r is  given  by  Bernoulli's 
equation  for  a compressible  fluid 


2 2 Y - 1 P 

Assuming  the  expansion  to  be 

— s constant  ,, 

pr 

and  that  H is  oonstant,  we  have  by  differentiating  (2)  and  (3)  that 

dV„  dVw  1 d_ 

Va  + Vw  — 2 + _ -E  = 0 
dr  W dr  p d,. 

and  using  (1) 

av„  dvw  v2 

VQ  — + V„— 2+—  =0  

a dr  w dp  r 

the  general  vortex  law  is 

vn  - c 

, Vw 

and  also  tan  as  — . . 

Va 

by  (5)  and  (6) 

Va  tan  arn*C  ..  
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By  differentiation 


n dVa  n_i  n 2 dn 

tan  a r — + Va  • tan  a n r + Va  • r seo  a — =0 


whence 


— = - i sin  2 a/—  ^ + £\ 
dr  2 \Va  dr  r/ 


(8) 


differentiating  (6) 


•w 


dVa 

2 

da 

_ 

— tan  a + 

Va  seo 

a — 

dr  " 

dj* 

a 

dr 

dVw 

CM 

5 

„ 2 

2 

11 

Va tan 

dr 

a + Va 

seo 

da 

d^ 


..  ..  (9) 


da 


* i * | - y . . , 

Substituting  for  — in  (9)  and  then  substituting  for  Vw and  — 

dr  dr  r 

in  (4)  finally  gives 


dV: 


Va  “ + (1  ~ n)  C r 


2 -(2n  + 1) 


= 0 


Integrating,  using  subscript  'o'  to  refer  to  conditions  at  the  inner 
radius  for  convenience 


2 2 
Va 


+ C 


2 (1  - n) 


Free  vortex  blading 

For  free  vortex  flow 
. . from  Equation  (10) 


n = 1 

V.  = V.  a oonstant 

(X  <Xq 


(10) 
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and  by  Equation  (7) 


J 0 

tan  a = — where  J = oonstant  = — 


2 „ C 

now  seo  a ■ 1 + — 5 

r d V i 
‘m  *a 


now  the  weighted  mean  radius 


2 2 ‘ 

= $ (H  + rc  ) 


o 2 2 

’ 2 „ 2 r0  tan  a0 

• “m  = 1 + —~2 2 

R + rQ* 


(11) 


Poroed  vortex  blading 


For  faroed  vortex  flow  n =>  - 1 


2 2 2 2 2. 
. , from  Equation  (10)  Va  = Va  - 2 C (r  - rQ  ) 


and 


tan  a = 


J\2  " 2 °2  (r2  " ro2) 


tan  am  = 


‘m 


r°  2/r2  _ 2\ 

o - 2 (rm  - r0  ) 


tan  a. 


seo  am  = 1 + 


2/_I° 


tan  ar 


, 2 2n 

(R  + rD  ) 

2 ~2  2 
- (R  - rc 


= 1 + 


(R2  ♦ rp2)  tan2  ap 
2\  r02  seo2  aQ  - R2  tan2  al 


= 1 + 


(R2  + rc2) 


2 (rQ2  ooseo2  aQ  - R2) 


..  (12) 
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